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ABSTRACT 





A seismic hazard analysis was developed for Guadalajara 
Metropolitan Area (ZMG), with a probabilistic approach, using the 
methodology developed by McGuire (2004). The ZMG is located on the 
central part of Jalisco State, in the Pacific Coast of Mexico. ZMG is comprised 
of six cities: Guadalajara, Zapopan, San Pedro Tlaquepaque, Tonalá, 
Tlajomulco de Zúñiga and El Salto. The seismic hazard settings were defined, 
using Uniform Hazard Spectra (UHS) and seismic parameters. Three Ground 
Motion Prediction Equations (GMPEs) were identified and a disaggregation 
study is presented. Geotechnical information was characterized for some 
ZMG sites; and hence site effects estimations were also investigated. UHS 
were established for both firm and soft soils for ZMG. Finally, the computed 
seismic parameters were compared to those for the seismic design guidelines 
for Guadalajara and Zapopan cities, and we found out that our UHS captured 
better the seismic settings of these cities of the Mexican Pacific Coast. 


Keywords: Seismic hazard, Uniform Hazard Spectra, site effects, seismic 
sources, GMPEs, recurrence relations, deaggregation, transfer function. 


ANÁLISIS DE PELIGRO SÍSMICO Y EFECTOS DE SITIO 
PARA LA ZONA METROPOLITANA DE GUADALAJARA 


RESUMEN 


En este trabajo se desarrolló un análisis de peligro sísmico para la 
Zona Metropolitana de Guadalajara (ZMG), con un enfoque probabilista, 
mediante la metodología desarrollada por McGuire (2004). La ZMG es 
localizada en la parte central del Estado de Jalisco en la Costa del Pacifico 
de México. La ZMG está conformada por seis municipios, Guadalajara, 
Zapopan, San Pedro Tlaquepaque, Tonalá, Tlajomulco de Zúñiga y El Salto. 
El peligro sísmico en la metodología probabilista se define con Espectros de 
Peligro Uniforme (EPU) y considera todas las fuentes sísmicas que afectan a 
la zona de estudio, se identificaron tres ecuaciones de movimiento y se 
presenta un estudio de desagregación. Se incluyen efectos de sitio a partir de 
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algunos estudios geofísicos de sitios ubicados en la ZMG. Se determinaron 
EPUs para terreno firme y para terreno blando de la ZMG. Los parámetros 
sísmicos de este trabajo, se compararon con la normatividad de diseño por 
sismo vigente de los municipios de Guadalajara y Zapopan, y encontramos 
mejores estimaciones sísmicas que las desarrolladas en estas normativas. 


Palabras clave: Peligro sísmico, Espectros de Peligro Uniforme, efectos de 
sitio, fuentes sísmicas, relaciones de atenuación, relaciones de recurrencia, 
desagregación, función de transferencia. 


1 INTRODUCCIÓN 


El peligro sísmico es una medida de la frecuencia de la ocurrencia de sismos 
con cierta intensidad; se cuantifica con la frecuencia con que excede un cierto nivel 
de aceleración o desplazamiento espectral y depende de la tectónica del territorio y 
de la sismicidad de las fuentes que intervienen en una región (Pérez et al., 2005). 
Un análisis de peligro sísmico envuelve una estimación cuantitativa de sismos 
ocurridos en un sitio en particular. Se obtiene mediante enfoques deterministas y 
probabilistas. Para este estudio se utilizó un enfoque probabilista con la metodología 
de Cornell (1968), desarrollada por McGuire (2004), la cual provee un marco de 
referencia donde las incertidumbres pueden ser identificadas, cuantificadas y 
combinadas. 


La metodología de Cornell (1968) utiliza dos modelos para la determinación 
del criterio del sismo de diseño a partir de un análisis probabilista: 


1. Un modelo de sismicidad, que describe la distribución geográfica de las 
fuentes de un evento y una distribución de magnitudes; que se expresa 
mediante relaciones de atenuación del movimiento fuerte del suelo. 


2. Un modelo de atenuación o GMPE que describe los efectos sobre un 
sitio dado mediante funciones que dependen de la magnitud, distancia 
de fuente a sitio, tipo de suelo, entre otras variables (Inca, 2013). 


La ecuación que se utiliza para evaluar el peligro sísmico de acuerdo a la 
metodología de Cornell (1968) es: 


N N 
MYX N 020=) y [| Piy>YI mre ham frio fate) dm ar de 
El El (1) 


Donde A(y > Y) es la tasa anual de excedencia del nivel del movimiento Y, 
debido a la ocurrencia de sismos en las N fuentes sísmicas que es igual a la suma 
de las tasas anuales de excedencia en cada zona fuente 2i(y > Y) las mismas que 
presentan una tasa anual de sismos vi; el término Pily > Y |m, r, e] da la probabilidad 
de excedencia condicionada al trio de variables m, r, e que representan la magnitud, 
distancia y epsilón; fMi (m)fRi (r) fei(e) son las funciones de densidad de 
probabilidad de la magnitud, distancia y epsilón (Aguiar y Rivas, 2018). 
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2 POTENCIAL SÍSMICO PARA LA ZONA METROPOLITANA DE 
GUADALAJARA 


En la evaluación del peligro sísmico es importante definir, los factores que 
influyen en la sismicidad de la región seleccionada, como aspectos geológicos y 
geotécnicos, y sobre todo se debe de estudiar el potencial sísmico de las posibles 
fuentes que afecten a la zona de estudio. 


Los grandes temblores a lo largo de la costa del pacífico en el Estado de 
Jalisco, son causados por la subducción de las placas oceánicas de Cocos y de 
Rivera, bajo la placa de Norteamérica y por ello son conocidos como sismos de 
subducción (ver figura 1). La placa de Rivera, se desplaza bajo la placa 
Norteamericana en el estado de Jalisco con una velocidad relativa de 2.5 cm/año 
frente a la costa de Manzanillo. La frontera entre las placas de Rivera y de 
Norteamérica se estima que interseca la costa de México cerca de Manzanillo 
(Ordaz, 2007). 


Placa del 
Pacífico 





Figura 1. Placas Tectónicas que interactúan en el Territorio Mexicano (SSN, 2016) 


En los últimos 160 años han ocurrido sismos muy importantes en el Bloque 
de Jalisco (ver figura 2), que han afectado a los estados de Jalisco, Colima y 
Michoacán con magnitudes representativas desde 5.1 hasta 8.2, desde el año 1567 
hasta el 2003, como el sismo del 3 de Junio de 1932, en la costa de Jalisco, donde 
se registró el sismo de mayor magnitud del siglo XX, que ha ocurrido en México, 
con una magnitud de 8.2 y réplicas de 7.8 y 6.9 los días 18 y 22 de junio de 1932, 
respectivamente (Nuñez-Cornú et al., 2004). 


Las principales estructuras tectónicas continentales más relevantes y de 
mayor impacto que se observan en el Bloque de Jalisco son: 
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Graben de Tepic-Zacoalco. Estructura con una orientación NW, con un 
corredor aproximado de 50 km por 250 km, es uno de los brazos del sistema 
triple de rift en el oeste de México. Junto con el rift de Colima y la Faja 
Volcánica Transmexicana, que limita el bloque de Jalisco (Martínez, 2011). 
El graben de Tepic-Zacoalco, generó sismos del 27 de diciembre de 1568 
con M = 7.5?, y el del 11 de febrero de 1875 con M = 7.5, pero una 
intensidad bastante grande de VIII (Lazcano, 2001). 


Graben de Chapala. Esta estructura contiene el lago natural de Chapala a 
1500 m de altitud, de 115 km de largo y de 10 a 25 km de ancho. El graben 
de Chapala se desarrolló en una falla transcurrente lateral izquierda dentro 
de la zona regional de falla Chapala-Tula E-W (Martínez, 2011). A éste 
graben se le asocia el sismo del 2 de octubre de 1847 con M=7.1 (Lazcano, 
2001). 


Graben de Colima. Estructura con anchura variable entre 20 y 60 km. El 
graben de Colima se divide en dos sectores, el graben norte, y el central. El 
graben está definido por fallas normales paralelas a su dirección, intersecta 
dos rifts: el de Zacoalco, al noreste, y al noroeste el rift de Citlala, para 
formar una unión triple, a unos 50 km al sur-suroeste de Guadalajara 
(Martínez, 2011). Es el graben con mayor actividad de los tres, entre los 
sismos más destacados por su magnitud y daños se destacan, los sismos 
de 1612 y 1785 con magnitudes menores a 6.0 y el del 22 de enero del 2003 
con Mw = 7.4 (Lazcano, 2001). 








Rvera 
SL 
gr! Raca del Pacífico 














Figura 2. Bloque de Jalisco (Núñez-Cornú, 2011) 
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2.1 Fuentes sísmicas tipo interplaca en la zona de subducción de 
México utilizadas 


Nishenko y Singh (1987) estudiaron los intervalos de recurrencia de sismos 
interplaca a lo largo de la zona de subducción de México (ver figura 3), de los trece 
segmentos o gaps característicos del país, seis se utilizaron en este estudio ya que 
se análizaron las fuentes interplaca que se encuentran dentro de un radio de 
influencia de 500 km, éste caso de estudio es explicado de manera detallada en el 
apartado 3.1. Los seis gaps utilizados se describen brevemente: 


1. Michoacán. Este segmento se encuentra en una longitud de 101.5W a 
103W. Coincide con la intersección batimétrica de la zona de fractura de 
Orozco con la parte media de la trinchera americana. La ocurrencia del 
sismo del 19 de septiembre de 1985, Ms = 8.1, indica que este gap tiene un 
modelo de grandes movimientos poco frecuentes. Se identificó un evento 
previo, el 7 de junio de 1911 con Ms = 7.9. El tiempo de recurrencia 
estimado para sismos tan grandes en esta región es de aproximadamente 
75 años. 


2. Colima l. Este segmento se encuentra en la longitud de 103 W a 103.7 W. 
Sismos previos asociados a este segmento son los del 15 de abril de 1941 
con Ms = 7.9 y el del 30 de enero de 1973 con Ms = 7.5. 


3. Colima Il. Este segmento se encuentra entre las longitudes de 103.7 "W a 
104.5"W. El gap de Colima es definido a base de las réplicas del sismo de 
Jalisco del 3 de junio de 1932 (Singh et al., 1985) y cubre una zona de 
aproximadamente 60 km de longitud entre la zona del sismo de 1932 y el 
de 1973. 


4. Jalisco. Entre las longitudes 104.3"W a 105.7 W. Los grandes sismos 
previos ocurridos en el estado de Jalisco fueron: el del 25 de marzo de 1806 
con Ms = 7.5; el del 31 de mayo de 1818 con Ms = 7.7; el del 20 de enero y 
16 de mayo de 1900 con Ms = 7.6 y 7.1 respectivamente y los del 3 y 18 de 
junio de 1932 con Ms = 8.1 y 7.8. El evento del 3 de junio de 1932 es uno 
de los eventos más grandes ocurridos en México hasta la fecha. Este sismo 
fue causado por una ruptura de la Placa de Rivera. 


5. Guerrero Central. Se encuentra entre las longitudes 100°W a 101 “W. Esta 
zona tiene la mayor tasa de recurrencia a lo largo de la zona de subducción. 
Sismos previos ocurridos en este segmento de Ms 2 7.5 fueron el del 7 de 
abril 18457, Ms = 7.9; 24 de diciembre de 1899, Ms = 7.7; 26 de marzo de 
1908, Ms = 7.8; 30 de julio de 1909, Ms = 7.5 y el del 16 de diciembre de 
1911, Ms = 7.8. 


6. Petatlán. Este segmento se encuentra entre las longitudes 101 W a 
101.8W. Sismos ocurridos en esta región con Ms 2 7.5 se incluyen los del 
22 de febrero de 1943, con Ms = 7.7; 14 de marzo de 1979, Ms = 7.9. En 
este segmento existe una gran variedad en los modos de ruptura. 
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Google earth 





Figura 3. Localización de los segmentos de ruptura de subducción tipo interplaca 
de la costa de México (Nishenko y Singh, 1987) 


Nishenko y Singh (1987), usaron los resultados de Nishenko y Buland 
(1987), quienes examinaron los intervalos de recurrencia para temblores 
característicos usando una función normalizada, T/Tave. Tave es el promedio del 
intervalo de recurrencia observado para un segmento de falla específico, y T es un 
intervalo de recurrencia individual. Basado en un conjunto de datos globales. 
Nishenko y Buland (1987) encontraron que la distribución logarítmica normal, 
proporciona un ajuste adecuado a la recurrencia de datos de tiempo para los 
segmentos de subducción, a diferencia de la distribución de Gauss o Weibull 
(Rikitake, 1976; Lindh, 1983; Sykes y Nishenko, 1984; y Nishenko, 1985). 


La función de densidad de probabilidad, f (T/T), definida por Nishenko y 
Buland (1987) se define por la siguiente ecuación: 


f(T/T)= —— et (T/T)-4D)?/20D? 


donde: 


T=es un intervalo de recurrencia individual 
oD=desviación estándar=0.215 
T=In(Taye) +uD 


Tave= promedio del intervalo de recurrencia observado para un segmento de falla específico 


uD = media = -0.0099 


El período de recurrencia esperado, Tex, es un número basado en 
observaciones de T para un segmento de falla específico, éste es definido como: 


= 2 
Tox = Terb 0D*/2 (3) 
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Siguiendo la terminología sugerida por Wallace et al. (1984), la estimación 
del tiempo de recurrencia anterior también se puede presentar en términos de una 
ventana de tiempo de pronóstico o predicción. La ventana utilizada en este estudio 
se define como el intervalo de tiempo. 


La ventana de tiempo usada en este estudio es definida por el intervalo de período: 


[Torea -n , Tprea+n ] 
(4) 


donde: 


Tpred= to+Texp 

Tprea = es el año del evento característico esperado obtenido del año del último sismo 
característico, más el período de recurrencia esperado 

to = es el año del último sismo característico 

n = 90% de intervalo de confianza para la estimación del Texp, donde: 


n=1.645(var[Toyp ]) 
var[Toxp ]=(5/ T?’+0D2)/T exp ? 


Los períodos de recurrencia aproximados por las ecuaciones (2), (3) y (4) aplicados 
a los segmentos de falla o fuentes de tipo interplaca se muestran en la tabla 1. 


Tabla 1. Período de recurrencia lognormal obtenida por Nishenko y Singh (1987), para los 
segmentos de falla interplaca 





Longitud | Ancho | to era n 
Segmento Ms 


km km año año | años 











1899 | 7.7 | 1960 39 





1908 | 7.8 | 1976 44 























Guerrero Central 139 80 
1909 | 7.5 | 1959 32 
1911 | 7.8 | 1979 44 
Petatlán 102 80 1979 | 7.6 | 2015 18 
Michoacán 133 80 1985 | 8.1 | 2015 37 
Colima 115 80 1973 | 7.5 | 2005 16 
Cas cal 2054 61 
e 180 80 | 1932 | 82 
Jalisco 


2010 51 
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Fuentes sísmicas tipo área o región utilizadas 


El territorio de México se subdividió, de manera global, en unidades 


sismogenéticas con el propósito de estimar el riesgo, basándose en los datos de los 
catálogos de registros tomados por instrumentos de medición e históricos (crónicas 
y escritos). El territorio de México se subdividió en 19 regiones principales (Zúñiga 
et al., 1997). Para nuestro estudio al igual que los segmentos interplaca, se 
analizaron las fuentes tipo área o región que se encuentran dentro de un radio de 
influencia de peligro sísmico de 500 km. Se tomaron en cuenta 11 fuentes tipo área 
o región, adicionales a los segmentos de falla explicados en el apartado 2.1 (ver 
figuras 4 y 5) de las cuales 10 son para sismos someros y 1 para sismos profundos. 
Las 11 regiones utilizadas se describen de forma breve: 


1. 


Región SUB1. Zona de eventos interplaca someros del tipo de subducción 
(profundidad < 40 km). Zona de transición de la convergencia entre las 
placas Rivera y Norteamericana (NOAM) a la convergencia Cocos - NOAM. 


Región SUB2. Zona de eventos interplaca someros del tipo de subducción. 
Convergencia Cocos - NOAM. 


Región RIV1. Zona de eventos interplaca someros de fallamiento normal 
principalmente, interface Pacífico — Rivera. 


Región RIV2. Zona de eventos interplaca someros de fallamiento de rumbo 
principalmente, interface Pacífico - Rivera. 


Región RIV3. Zona de eventos interplaca de subducción somera. Interface 
Rivera - NOAM. 


Región BB. Zona de eventos intraplaca (NOAM) someras, provincia de la 
Cuenca de Burgos (Burgos Basin). 


Región MVB. Zona de eventos intraplaca (NOAM) someros (h < 40km). 


Región NAM. Esta zona es de eventos intraplaca (NOAM) someros. Esta 
región comprende la actividad somera (profundidad menor a 40 km). 


Región BAR. Zona de eventos intraplaca (NOAM) someros, provincias de 
Cuencas y Sierras - Fisura del Río Bravo. 


. Región BC1. Esta zona comprende la actividad asociada a los eventos 


someros intraplaca en el área de Baja California. 


. Región IN1. Zona de eventos intraplaca (Cocos) de profundidad intermedia 


(40 km <h < 120 km). Relacionados a la interface Cocos - NOAM. 


Zúñiga et al. (1997), propone idealizar el proceso de ocurrencia de 


temblores como uno de poisson. En éste proceso los tiempos entre eventos se 
encuentran de forma exponencial y puede demostrarse que su distribución no se 
altera al aumentar el tiempo transcurrido sin un temblor. 
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En este proceso la densidad de probabilidades del tiempo entre eventos con 
magnitud mayor o igual a M es exponencial y está dada por: 


P(t) = Me * (5) 


Donde t es el tiempo entre eventos con magnitud mayor que M y A(M) es la 
tasa de excedencia de la magnitud M, definida como el valor esperado del número 
de temblores con magnitud mayor o igual M, por unidad de tiempo. Si se conoce 
A(M), la distribución de los tiempos entre eventos está totalmente determinada. A(M) 
es una función que decrece con la magnitud, y se le ha asignado la siguiente forma 
funcional (Cornell y Vanmarcke, 1969): 


oPM_g-BM, 


MM) = Mo e-BMo-e-BMy (6) 


donde: 


Ao, B y Mu = son parámetros que definen la tasa de excedencia de cada una de las 
fuentes sísmicas 

Mo = es la magnitud por encima de la cual el catálogo esta completo, para todas se 
considera 4.5 


Éste modelo teórico fue adoptado para los temblores costeros con M<7, y 
para el resto de los temblores. Para los temblores superficiales, el proceso real de 
ocurrencia fue dividido artificialmente en dos, teniendo como límite la magnitud M = 
7. Esto implica que, ésta es la magnitud máxima posible en la parte del proceso 
considerada poissoniana (Zúñiga et al., 1997). Los parámetros de recurrencia 
obtenidos para las regiones tectónicas tipo región o área usadas en este estudio se 
pueden observar en la tabla 2. 


Tabla 2. Parámetros de períodos de recurrencia de las regiones o zonas tectónicas tipo 
área consideradas en el estudio (Zúñiga et al., 1997) 
































E(Ao) 
Zona E(B) c(B) E(Mu) Mo 
1/año 
SUB 1 0.041 2.42 0.282 7.2 4.5 
SUB 2 3.484 2.210 0.089 7.2 4.5 
RIV 1 0.595 2.390 0.195 7.2 4.5 
RIV 2 4.512 2.390 0.079 7.2 4.5 
RIV 3 0.191 2.660 0.248 7.2 4.5 
BB 0.080 2.410 0.286 5.4 4.5 
MVB 0.409 2.480 0.210 7.2 4.5 
NAM 1.756 3.160 0.076 7.4 4.5 


























98 Angélica Buenrostro, Alonso Gómez y Hugón Juárez 

















E(Ao) 
Zona E(B) c(B) E(Mu) Mo 
1/año 
BAR 0.083 2.370 0.280 5.0 4.5 
BC 1 0.963 2.490 0.199 6.0 4.5 
IN 1 6.924 2.410 0.061 7.9 45 





























Figura 4. Fuentes tipo área para sismos someros interplaca e intraplaca usados 
para los casos estudiados (Zúñiga et al., 1997) 





Figura 5. Regiones sismotectónicas de subducción tipo intraplaca para sismos 
profundos (Zúñiga et al., 1997) 
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2.3 Modelos de recurrencia 


Las relaciones de recurrencia para las fuentes sísmicas utilizadas, fueron 
las propuestas por Nishenko y Singh (1987) para fuentes tipo interplaca. Y las de 
Zúñiga et al. (1997) para eventos de tipo área o región, desarrolladas en los 
apartados 2.1 y 2.2 del presente capítulo. Ambas propuestas se incluyen 
actualizadas en el programa de cómputo empleado para el análisis de peligro 
sísmico, EZ-FRISKv7.65'M, el cual arroja estas relaciones como un modelo 
Gutenberg — Richter Truncado. 


En la figura 6, se muestra la relación de la recurrencia anual de las fuentes 
sísmicas usadas, con la magnitud de momento (Mw) que pueden tener los eventos 
causados por estas fallas, la figura nos indica un arreglo del modelo magnitud- 
recurrencia, Gutenberg-Richter Truncado. 
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Figura 6. Relación de Recurrencia de las fuentes sísmicas utilizadas para un 
r=500 km con un modelo Gutenberg-Richter Truncado (Buenrostro, 2017) 


2.4 Modelos de atenuación (GMPE;) 


Una vez determinada la tasa de actividad de cada una de las fuentes 
sísmicas, es necesario evaluar los efectos que, en términos de intensidad sísmica, 
produce cada una de ellas en un sitio de interés. Para ello se requiere saber que 
intensidad se presentaría en el sitio, si en la i-ésima fuente ocurriera un temblor con 
magnitud dada. A las ecuaciones que relacionan magnitud, posición relativa de la 
fuente-sitio e intensidad, se les conoce como leyes de atenuación (Ordaz et al., 
2007). Para este estudio se seleccionaron 3 leyes de atenuación de acuerdo a las 
fuentes caracterizadas para nuestro análisis de peligro sísmico, se involucraron 
modelos en los cuales se utilizaron catálogos de sismos de México, que 
representaran principalmente una respuesta espectral adecuada de los sismos de 
costa o subducción. Ya que el análisis de peligro sísmico no se ésta realizando para 
una edificación en particular, los modelos proporcionan suficiente información para 
encontrar espectros de respuesta elástica para un rango de período de interés 
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(Aguiar y Rivas, 2018). Los tres modelos utilizados son Youngs et al. (1997), García 
et al. (2005) y Arroyo et al. (2010). 


2.4.1 Youngs et al. (1997) 


Este modelo fue desarrollado por RR. Youngs, S.J. Chiou, W.J.Silva y J.R. 
Humphrey. Se presenta una relación de atenuación para sismos de subducción, y 
sismos provocados por fallas de interplaca e intraplaca, con magnitudes de 
momento mayores a 5.0 y distancias de 10 a 500 km. Esta ley es empleada en el 
análisis ya que usa una base de datos de registros internacionales y dentro de ella 
se encuentran 23 registros de sismos fuertes en la zona de subducción de México. 
Las relaciones de atenuación fueron desarrolladas para aceleraciones pico (PGA) 
de zonas de subducción que son evaluadas mediante análisis de regresiones 
aleatorias de datos empíricos. Las constantes son determinadas por un análisis que 
incluye coeficientes de las relaciones de atenuación seleccionados para grupos de 
referencia y parámetros que representan perturbaciones de otros 5 grupos de 
referencias, C1 a C5 (Youngs et al., 1997). El modelo está limitado a analizar hasta 
4 segundos la respuesta espectral elástica, por lo cual, para edificaciones de 
períodos altos no es conveniente utilizarla, sin embargo, como se explicó en el 
párrafo anterior, éste modelo se adopta para éste estudio, ya que no se analiza el 
peligro sísmico de un edificio en específico sino se realiza un estudio general para 
la ciudad. 


La ecuación que define esta ley para roca está dada por: 


In SA(T)=0.2418+1.414 M+C,(T)+C2(T)*(10-M)? +C3(T)* In(rrup+1.7818 
20.55 M) +0.00607 H+0.3836 ZT (7) 


Desviación estándar=C¿+Cs5 M (8) 
donde: 


SA (T) = aceleración espectral normalizada en g, calculada para el período T en 
estudio 

T = período en segundos 

C1(T),C2(T),C3(T),C4(T),C5(T)= dependientes del período T, obtenidos por medio 
de regresión de datos dados por Youngs et al. (1997) 

M = magnitud de momento 

H = profundidad al foco en km 

ZT = tipo de falla; usar O para falla interplaca y 1 para falla intraplaca 

rrup = distancia más cercana a la ruptura en km o distancia hipocentral. El intervalo 
de distancias es de 10 a 500 km 


2.4.2 García et al. (2005) 


Por D. García, S.K. Singh, M. Herraíz, M. Ordaz y J.F. Pacheco. Este 
modelo fue desarrollado para sismos inslab o intraplaca de las costas de México 
(falla normal, profundidad intermedia), para sitios sobre roca firme, el número de 
eventos sísmicos que se utilizaron fueron 16 a distancias locales y regionales de 
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hasta 400 km, con momentos de magnitud (Mw) entre 5.2 y 7.4 y profundidades 
focales (H) entre 35 y 138 km (García et al., 2005). 
La ecuación que define esta ley para roca se define como: 


log SA(M=C+(M+C2(T) My +C3(T) R-C4(T) log(R) +C5(T) H+0 (9) 
donde: 


SA(T) = aceleración espectral horizontal para 5% de amortiguamiento en cm/s? 

T = período en segundos 

Mw = magnitud de momento 

H = profundidad al foco en kilómetros 

C¿M*C1M+C(M+CM+Cs(TM) = constantes definidas por un análisis de 
regresión al conjunto de datos dependientes de cada período en García et al. (2005). 
R = es la distancia promedio de la distancia a la superficie de falla en km, definida 
por R = (R2c/d +A2)1/2, donde Re/d, es la distancia más cercana a la superficie de 
falla para eventos con Mw > 6.5, o la distancia hipocentral para el resto y A = 
0.0075x100.507 Mw, el cual es un término de saturación de la fuente definido por 
Atkinson y Boore (2003) 

O = desviación estándar 


2.43 Arroyo et al. (2010) 


Este modelo se desarrolló por D. Arroyo, D. García, M. Ordaz, M. A. Mora y 
S. K. Singh. Se derivaron relaciones de atenuación para componentes horizontales 
de pseudoaceleraciones de espectros de respuesta para sismos interplaca en 
México en roca (Clase B en la clasificación NEHRP) en la región de la costa. El 
modelo resultante tiene una precisión similar a los modelos derivados de bases de 
datos regionales y mundiales de zonas de subducción. Los resultados de este 
modelo revelan que la relación de atenuación en México tiende a ser más fuerte que 
las relaciones de atenuación mundiales para grandes eventos; y predicen éstas 
relaciones para grandes sismos a grandes distancias de la fuente a diferencia de 
cómo se podrían considerar para modelos globales. 


Para un modelo circular de fuente finita en zonas de subducción la ecuación 


se define como: 
E; (a4(T) R)-E4 (aT) VR? +ro?) 


In SA (M= a¿(D+az,CT) My+as (In ro? 
(10) 


donde: 


SA(T) = aceleración espectral horizontal para 5% de amortiguamiento calculada 
para T 

ai(T) = son coeficientes determinados atreves de un análisis de regresión 

R = distancia más cercana a la superficie de falla en km 

ro = radio equivalente a una falla circular y se define como: ro? = 1.4447 x 10-5 e 
2.3026 Mw 

Mw= magnitud de momento 

Ey = integral exponencial, obtenida mediante métodos numéricos dada por: 


vo pt 
E= f ot 
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3 ANÁLISIS DE PELIGRO SÍSMICO 


3.1 Caso de estudio para diferentes radios de influencia 


Se realizó un análisis para determinar el radio de influencia de peligro 
adecuado para la zona de estudio. Se obtuvieron Espectros de Peligro Uniforme 
(EPU), calculados para períodos de retorno de 20 años, 45 años, 100 años, 475 
años, 975 años y 2475 años, que representan probabilidades de excedencia 
respectivamente, del 91.8%, 67.1%, 39.3%, 10%, 5%, 2%, de ocurrir en 50 años 
(vida útil de una edificación). Se localizaron los centros históricos de cada uno de 
los 6 municipios de la ZMG, (ver figura 7 (a)) y se trazaron radios de influencia de 
200 km, 300 km y 500 km para cada uno de estos puntos de localización. Los 
círculos de influencia se definen para observar y estudiar las variantes de las 
diferentes fallas que se encuentran en estos radios, de manera que se pueda discutir 
la influencia de los segmentos y áreas de falla sobre los Espectros de Peligro 
Uniforme (Inca, 2013). En la tabla 3, se observan las fuentes sismogénicas y las 
características utilizadas para cada caso de estudio. 


Tabla 3. Datos utilizados para el análisis de peligro sísmico de los casos de 
estudio para diferentes radios de influencia (Buenrostro, 2017) 



























































Tipo Fuentes sismogénicas Modelos de 
Nombre | Municipio Latitud Longitud de Fuentes 7 A Pes 
Suelo Falla Fuentes Área OMEGN 
IN1-C, IN1-B, IN- YOU-1997, 
D Guadalajara | 20.6766° | -103.3467° Firme SINT A, SUB1, GA-2005, 
y SUB2, BB, MVB ARR-2010 
IN1-C, IN1-B, IN- 
; y YOU-1997 
GDL- : A 3 F JAL, COL, A, SUBA, f s 
300 Guadalajara | 20.6766 -103.3467 Firme MI-CO, MI RIV3 SUB2, BB, R 
MVB, NAM E 
GDL- JAL CoL aya | YOU-1997, 
Guadalajara | 20.6766* | -103.3467* | Firme MI-CO, MI, i ñ , GA-2005, 
300 PE, GRO-C Ped ARR-2010 
i È NAM,RIV1,RIV2 A 
IN1-C, IN1-B, IN- YOU-1997, 
ea Zapopan 20.6866* | -103.3508* Firme Aar A, SUB1, GA-2005, 
É SUB2, BB, MVB ARR-2010 
IN1-C, IN1-B, IN- 
i ? YOU-1997 
ZAP- o E j JAL, COL, A, SUB1, RIV3, a : 
300 Zapopan 20.6866 -103.3508 Firme MI-CO, MI SUB2, BB, MVB, GA-2005, 
ARR-2010 
NAM 
Ta JAL CoL o ua | YOU-1997, 
Zapopan 20.6866° | -103.3508° | Firme MI-CO, MI, ; A GA-2005, 
200 PE, GRO-C SOBERE PNVD: ARR-2010 
+ 7 NAM,RIV1,RIV2 s 
IN1-C, IN1-B, IN- YOU-1997, 
eS Tlaquepaque | 20.6354* | -103.3053* Firme se AL A, SUBI, GA-2005, 
: SUB2, BB, MVB ARR-2010 
IN1-C, IN1-B, IN- 
A i YOU-1997 
TLAQ- $ 5 ; JAL, COL, A, SUBA, E i 
300 Tlaquepaque | 20.6354 -103.3053 Firme MI-CO, MI RIV3 SUB2, BB, a 
MVB, NAM Ñ 
IN1-C, IN1-B, IN- 
TLAQ- JAL, COL, A, SUBA, YOU-1997, 
500 Tlaquepaque | 20.6354* | -103.3053” | Firme MI-CO, MI, RIV3,SUB2, BB, GA-2005, 
PE, GRO-C MVB, ARR-2010 
NAM,RIV1,RIV2 
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Tipo Fuentes sismogénicas Modelos de 
Nombre Municipio Latitud Longitud de Fuentes y A me 
Suelo Fuentes Área tenuacion 
Falla 
IN1-C, IN1-B, IN- YOU-1997, 
as Tonalá 20.6473° | -103.2673° | Firme Aa A, SUBA, GA-2005, 
’ SUB2, BB, MVB ARR-2010 
IN1-C, IN1-B, IN- 
: l YOU-1997 
TON- i A š . JAL, COL, A, SUB1, RIV3 i 
300 Tonalá 20.6473 -103.2673 Firme MI-CO, MI SUB2, BB, MVB, GA-2005, 
NAM ARR-2010 
IN1-C, IN1-B, IN- 
TON- JAL, COL, A, SUB1, RIV3, YOU-1997, 
500 Tonalá 20.6473° | -103.2673° Firme MI-CO, MI, SUB2, BB, GA-2005, 
PE, GRO-C | MVB,NAM, RIV1, ARR-2010 
RIV2 
IN1-C, IN1-B, IN- YOU-1997, 
Le Tlajomulco 20.4737* | -103.4479* Firme a A, SUBI, GA-2005, 
i SUB2, BB, MVB ARR-2010 
IN1-C, IN1-B, IN- 
; , YOU-1997 
TLAJ- A 5 $ , JAL, COL, A, SUB1, RIV3, ` ? 
300 Tlajomulco 20.4737 -103.4479 Firme MI-CO, MI SUB2, BB, MVB, GA-2005, 
NAM ARR-2010 
IN1-C, IN1-B, IN- 
TLAJ- JAL, COL, A, SUB1, RIV3, YOU-1997, 
500 Tlajomulco 20.4737* -103.4479* Firme MI-CO, MI, SUB2, BB, GA-2005, 
PE, GRO-C | MVB,NAM, RIVA, ARR-2010 
RIV2 
IN1-C, IN1-B, IN- YOU-1997, 
po El Salto 20.5196° | -103.1784° | Firme CSA A, SUBA, GA-2005, 
dl SUB2, BB, MVB ARR-2010 
IN1-C, IN1-B, IN- 
i j YOU-1997 
SAL- i à ; JAL, COL, A, SUB1, RIV3, ; 
300 El Salto 20.5196 -103.1784 Firme MI-CO, MI SUB2, BB, MVB, GA-2005, 
NAM ARR-2010 
IN1-C, IN1-B, IN- 
SAL- JAL, COL, A, SUB1, RIV3, YOU-1997, 
n o | Firme | MI-CO, MI, SUB2, BB, GA-2005, 
500 El Salto 20.5196 -103.1784 PE, GRO-G | MVB,NAM, RIVA, ARR-2010 
RIV2 
































En las figuras 7 (b), 7 (c) y 7 (d) se muestran el trazo de estos radios, para 
el centro de Guadalajara, uno de los puntos seleccionados para este caso de 
estudio. Se obtuvieron para cada uno de los radios de influencia EPUs para los 
períodos de retorno mencionados en el párrafo anterior. 


En la figura 8 se muestra una comparativa de los tres radios de influencia 
para los puntos del centro de Guadalajara, GDL-200, GDL-300 y GDL-500. En esta 
comparación se puede observar que para lo radios de 200 km y 300 km los 
resultados se sobreponen, ya que una o algunas fallas que se encuentren a menos 
de 200 km puedan ser las que ocasionen un peligro mayor a la zona estudiada. En 
los EPUs obtenidos para los 500 km aumentan ligeramente las ordenadas 
espectrales, en los períodos de O a 0.5 s; ya que la contribución de las demás 
fuentes es muy pequeña. Se tomó la decisión de utilizar un radio final de 500 km, 
que abarcará un mayor número de fuentes, para el análisis de peligro sísmico de 
cada uno de los puntos de la malla de la microzonificación de la ZMG. 
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Figura 7. (a) Localización de los puntos para la selección de radios de influencia 
(b) Radios de influencia para segmentos de falla (c) Radios de influencia para 
áreas de sismos someros (d) Radios de influencia para áreas de sismos profundos 
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R-200 T=20 años 
R-200 T=45 años 
R-200 T=100 años 
R-200 T=475 años 
R-200 T=975 años 
R-200 T=2475 años 
R-300 T=20 años 
R-300 T=45 años 
R-300 T=100 años 
R-300 T=475 años 
R-300 T=975 años 
R-300 T=2475 años 
+ = R-500 T=20 años 
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Comparación de EPU GDL-200, GDL-300 y GDL-500 
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Aceleración Espectral SA (g) 





Figura 8. Espectros de peligro uniforme para tres radios en el punto “Guadalajara”, 
las curvas R=200 km y R=300 km (líneas sólidas) se superponen (Buenrostro, 
2017) 


3.2 Malla de interpolación (IDW) 


Se realizó una malla de puntos con el fin de obtener una microzonificación 
de la ZMG y calcular los EPUs a una distancia específica en cada una de estas 
intersecciones. Para elegir esta distancia se basó en la investigación realizada por 
Inca (2013), que mediante el método de Interpolación de Ponderación de Distancia 
Inversa (IDW), de Shepard (1968) concluyó, que se podía tener una buena 
aproximación de los EPUs desconocidos, con distancias menores o iguales a 15 y 
30 km debido a las cortas distancias entre ellas. 


Para este estudio se tomó una distancia arbitraria de 10 km debido a que la 
extensión territorial de la ZMG no es muy grande. Por lo que se decidió refinar esta 
malla de esa forma y en consecuencia obtener resultados con buenas 
aproximaciones para la interpolación. En la figura 9 (izquierda), se observa esta 
malla, donde se localizaron 62 puntos para poder trazar la cuadricula. Se 
identificaron los puntos que se encuentran dentro de los 6 municipios de la ZMG, 
(zona de color gris oscuro), quedando así, un total de 24 puntos para este estudio, 
como se observa en la figura 9 (derecha). 


A estos puntos se les asignó una clave para distinguirlos y se obtuvieron 
sus ubicaciones exactas para el análisis (latitud y longitud). Se tomó un radio de 
influencia de fuentes de fallas de 500 km. Para cada punto se le dio una relación de 
atenuación a cada falla seleccionada dentro del radio de influencia, como en los 
casos de estudio desglosados en la tabla 3. Y, por último, se consideró una 
velocidad de onda de corte del suelo a 30m (Vs30), mayor o igual a 760 m/s, para 
suponer un terreno firme o tipo roca, según los códigos americanos como el ASCE 
7-16, para cada uno de los puntos, ya que los modelos de atenuación usados están 
calculados para esas características de suelo. Posteriormente, se aplicarán efectos 
de sitio sobre los EPUs obtenidos mediante estudios geofísicos (en el capítulo 4 se 
explica de forma detallada). 
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Figura 9. Trazo de malla de interpolación con intersecciones a cada 10 km 
(izquierda) y malla de interpolación definitiva con 24 puntos (derecha) (Buenrostro, 
2017) 


En las figuras 10 y 11 se muestran los Espectros de Peligro Uniforme, 
obtenidos mediante el análisis de peligro sísmico, realizado para los puntos: ZMG- 
01 y ZMG-17. Como representación de los valores espectrales mayores y menores 
de los 24 puntos obtenidos en la malla de interpolación. 












1,00 
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0,00 0,50 1,00 1,50 2,00 2,50 3,00 
Período T (s) 


Figura 10. EPUs, obtenidos para el punto ZMG-01 (respuesta mayor), ubicado en 
el municipio de Tlajomulco (Buenrostro, 2017) 


Las aceleraciones máximas para todo los EPUs, se ubican en los períodos bajos 
(Os a 0.55), ya que se obtuvieron para terreno firme o roca, donde aproximadamente 
este tipo de suelo corresponde a períodos menores a 0.5 segundos. Se puede 
observar que el punto ZMG-01 (figura 10), es donde se tienen mayores 


Análisis de Peligro Sísmico y efectos de sitio... 107 


aceleraciones espectrales respecto a los 24 puntos calculados). Este punto se 
encuentra en Tlajomulco, y es el más cercano a las costas de Colima y de Jalisco. 
Por lo que, es muy probable que ese sea el motivo de las ordenadas de 
pseudoaceleración de mayor amplitud, y que se tenga un mayor peligro sísmico 
provocado por la subducción de la costa. En el EPU obtenido para el punto ZMG- 
17 en Zapopan (ver figura 11), las aceleraciones espectrales son menores que los 
demás puntos mostrados, ya que es el punto más alejado al Pacífico. 
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Figura 11. EPUs, obtenidos para el punto ZMG-17 (respuesta menor), 
correspondiente al municipio de Zapopan (Buenrostro, 2017) 
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Aceleración espectral (9) 


Figura 12. Curvas de peligro por fuente sísmica para el punto Z2MG-01, para un Ts 
= 0.3 s (Buenrostro, 2017) 


Se muestran las contribuciones que cada falla aporta (ver figura 12 y 13), al peligro 
sísmico de los puntos ZMG-01 y ZMG-17 para un período de 0.3 segundos (período 
en donde la SA es mayor). En ambos casos se aprecia que las fallas que aportan 
un mayor peligro, son las fallas interplaca del Gap de Colima y de Jalisco y la falla 


108 Angélica Buenrostro, Alonso Gómez y Hugón Juárez 


tipo área IN1-A (curvas con líneas punteadas en las gráficas); siendo los eventos 
de subducción los más importantes a considerar para este análisis. 
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Figura 13. Curvas de peligro por fuente sísmica para el punto Z2MG-17, para un Ts 
= 0.3 s (Buenrostro, 2017) 


3.3  Desagregación del peligro sísmico 


La desagregación del peligro sísmico, es una manera efectiva de identificar 
diferentes escenarios de eventos (expresados en términos de distancia — magnitud), 
que contribuyen al peligro sísmico de un sitio. El escenario de la contribución de los 
eventos al peligro sísmico de un sitio, C (M, R), puede ser obtenido mediante la 
determinación de la distancia a la fuente (r) desde el sitio, (Filiatrault et al., 1998). 
Este proceso muestra la contribución a la frecuencia anual de excedencia de 
acuerdo a la magnitud, distancia y parámetros de desviación estándar de las 
relaciones de atenuación. Si el valor de aceleración ocurre en el lugar de interés 
para la evaluación de la amenaza, ciertas magnitudes, distancias y valores de 
desviación contribuyen a dicha aceleración más que otros. El proceso de 
desagregación muestra estas contribuciones relativas. Es posible analizar la 
variación de tres parámetros magnitud, distancia y densidad de probabilidad 
(Kramer, 1996). 


Se realizó la desagregación del peligro sísmico para cada uno de los 24 
puntos de la malla, para un r = 500 km, para un período de retorno Tr = 475 años, 
ya que es el período de retorno asociado a un espectro de diseño básico para 
estructuras tipo B, de acuerdo a la normativa de los municipios de Guadalajara y 
Zapopan. Se consideraron períodos de vibración de 0.2 s, 0.75 s y 2.0 s. Se eligieron 
estos períodos por el número de niveles promedio que tienen aproximadamente la 
mayoría de las edificaciones que se encuentran en la ZMG (2 a 20 niveles). En las 
figuras 14 y 15 se muestran los resultados de la desagregación con los parámetros 
mencionados en el párrafo anterior, para los puntos ZMG-01 y ZMG-17. Se muestra 
la desagregación de estos dos puntos debido a que son los que tienen valores de 
aceleraciones espectrales, mayores y menores respecto a los demás puntos 
analizados. 
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Figura 14. Desagregación del peligro sísmico (M y R) para el punto ZMG-01, para 


períodos 0.2 s, 0.75 s y 2 s (Buenrostro, 2017) 


Se observa una magnitud promedio de 8.4 - 8.65. Las distancias promedio que 
arrojan para los tres puntos son diferentes por la localización en la que se 
encuentran. Por el tipo de magnitud que indica, y las distancias promedio, las 
fuentes sísmicas que están ocasionando un mayor peligro a la ZMG, son las de las 
costas de Jalisco y de Colima. Estos resultados sirvieron para corroborar lo que se 
había observado en las curvas de peligro por fuente sísmica. En todos los puntos 
de la malla se observaron resultados muy similares, por lo que estas características 
particulares producto de la desagregación, nos ayudaron, en la búsqueda de 


acelerogramas para los efectos de sitio aplicados en los EPUs (la metodología 


usada es explicada en el capítulo 4). 
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Figura 15. Desagregación del peligro sísmico (M y R) para el punto ZMG-17, para 
períodos 0.2 s, 0.75 s y 2 s (Buenrostro, 2017) 


4 ¿ESPECTROS DE PELIGRO UNIFORME CON EFECTOS 
LOCALES 


La desagregación del peligro sísmico obtenida, nos sirvió como parámetro 
para la búsqueda de acelerogramas, que estos a su vez, nos ayudarán a obtener la 
amplificación de los efectos locales, de las ondas de cuerpo del suelo de la zona 
estudiada. Se analizaron 14 registros de la ZMG, del sismo del 9 de octubre de 1995 
ocurrido en Manzanillo, Colima. Los registros fueron obtenidos de la Base Mexicana 
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de Sismos Fuertes (BMSF) volumen 2, la cual tiene los acelerogramas registrados 
en la República Mexicana desde el año, 1960 a 1999. También se encontraron y 
analizaron 11 registros en el Instituto de Ingeniería de la UNAM (I de | UNAM) para 
la ZMG, para los años del 2012 al 2015. 


Se procesaron todos los datos, para las tres componentes de cada uno de 
los 25 registros (vertical, horizontal N-S y horizontal E-O) y se corrigieron por línea 
base, con un filtro tipo Butterworth con una frecuencia baja de 0.10 Hz y una 
frecuencia alta de 25 Hz, para obtener espectros de respuesta con un 5% de 
amortiguamiento. 


Aceleración N90W 
Sismo del 09 de octubre de 1995, Cicej-Pozo 35 m 
(a) 
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Sismo del 09 de octubre de 1995, Cicej-Pozo 35 m 
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Figura 16. Acelerogramas elegidos, registrados en tres diferentes componentes: 
horizontal N-S (11.a), vertical (11.b) y horizontal E-O (11.c) 


Las características que se tomaron en cuenta para elegir un registro, y 
utilizarlo en la obtención de efectos locales del suelo, fueron las siguientes: 
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e  Acelerogramas registrados en alguna estación acelerográfica de la Zona 
Metropolitana de Guadalajara. 

e Registros de sismos con magnitud > 7. 

e  Acelerogramas de sismos ocurridos en la costa del Pacífico, con 
características sísmicas similares, a las del sismo del 3 de junio de 1932. 

e  Acelerogramas registrados en terreno firme o roca. 


Debido a que la instrumentación sísmica de la ciudad de Guadalajara es 
escasa, se tomó la decisión de utilizar el registro de la estación: CICEJ Pozo 35m, 
del sismo 9 de octubre de 1995, el cual reunió con las características deseadas. En 
la figura 16, se muestran los acelerogramas del registro elegido en las tres 
componentes medidas. 


4.1 Funciones de transferencia 


Una función de transferencia es la relación de amplificación que hay entre 
la base (terreno firme o roca) y la superficie libre. Para tener como resultado 
funciones de transferencia de nuestros suelos, se requieren ecuaciones de onda, 
que puedan obtener la amplificación entre 2 estratos, sus deformaciones y 
aceleraciones (Gómez-Bernal et al., 2007). 


Para esta investigación, se desea obtener la respuesta del suelo 
suponiendo propagación vertical de ondas de cortante, para este fin se utilizó el 
programa de cómputo, SHAKE-91_input (Idriss & Sun, 1992), versión 2013, que 
utiliza el método lineal equivalente. Este método, puede tener una aproximación 
acertada de los efectos no lineales del suelo en la respuesta, llevando a cabo una 
secuencia de análisis lineales, que consideran valores constantes del módulo de 
rigidez de la relación de amortiguamiento. Para obtener la amplificación del suelo, 
se modeló cada estrato de cada uno de los 18 perfiles de suelo del poniente de la 
ZMG, los cuales fueron proporcionados por Lazcano (2012). Se modeló, la 
velocidad de onda de corte (Vs), peso volumétrico (y) y amortiguamiento mínimo, 
se le asignó una curva de degradación de rigidez y amortiguamiento para cada 
estrato. Las curvas usadas fueron las de Seed & Idriss (1970), para arenas 
promedio. Los 18 perfiles de suelo y sus características geofísicas obtenidas con un 
estudio de REMI (Refracción de Microtremores) se muestran en la tabla 4. 


Tabla 4. Perfiles de suelo (Lazcano, 2012 y 2016) 










































































Número Ubicación Municipio Estratos 1 2 3 4 5 6 
Profundidad (m) | 0-3 3-6 6-9 9-22 |22-32| >32 
Vs (m/s, 152 186 396 262 443 880 
1 Catedral de Guadalajara (ms) 
Guadalajara y (kN/m3) 12.60 | 14.10 15.00 14.10 | 16.60 | 17.20 
T (s) 30m 0.47 
, Profundidad (m) | 0-4 4-9 9-25 >25 
agua W. /: 180 123 309 1300 
ibli s (m/s 
2 Biblioteca Guadalajara (ms) 
TE i y (kN/m3) 14.10 | 12.40 14.80 17.40 
(Agua Azul) Ts (s) 30m 0.45 
Profundidad (m) | 0-1.5 | 1.5-10 | 10-32 > 32 
i ivi Vs (m/s 143 298 363 1260 
3 Registro Civil Guadalajara (m/s) 
No.1 y (kN/m3) 12.50 | 14.80 | 14.80 17.30 
Ts (s) 30m 0.41 
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Número Ubicación Municipio Estratos 1 2 3 4 5 
Profundidad (m) | 0 - 4.5 | 45-11 | 11-50 > 50 
Fco. Javier ) Vs (m/s) 154 369 382 850 
4 Gamboa y Lerdo | Guadalajara 
de Tejada y (kN/m3) 12.6 | 147 | 1470 | 17.20 
Ts (s) 30m 0.58 
Profundidad (m) | 0-5 5-13 13-53 | 53-54 > 54 
Hotel Riu Plaza, ] Vs (m/s) 180 290 450 235 950 
5 López Mateos y | Guadalajara 
Lázaro Cárdenas Y (kN/m3) 1400 | 148 | 16:10 | 1480 | 17.20 
Ts (s) 30m 0.60 
Profundidad (m) | 0-5 5-11 11-57 > 57 
6 Centro i z Vs (m/s) 165 290 575 1000 
comercial, apopan 
La Gran Plaza y (kN/m3) 132 | 148 | 16.30 | 17.30 
Ts (s) 30m 0.52 
Profundidad (m) | 0-3 3-12 12-44 > 44 
i Vs (m/s, 167 217 596 1200 
A Patria y Zapopan (m/s) 
Guadalupe y (KkN/m3) 13.2 14.5 16.40 17.30 
Ts (s) 30m 0.45 
A A Profundidad (m) | 0-6 6-11 11 -24 24 - 62 >62 
vii Vs (m/s) 175 260 420 610 860 
Panamericana s (m/s, 
8 Zapopan 
Campus Gdl, Cd. de y (kN/m3) 13.4 | 146 16.10 | 16.50 | 17.00 
Granja Ts (s) 30m 0.63 
Profundidad (m) | 0-3 3-11 11-25 | 25-85 > 62 
ó López Mateos y | zapopan Vs (m/s) 185 260 445 450 800 
Mariano Otero Y (kN/m3) 13.40 14.60 16.20 16.20 17.00 
Ts (s) 30m 0.80 
Profundidad (m) | 0-3 3-9 9 - 31 > 31 
Colomos y ñ Vs (m/s) 188 253 399 1400 
10 Manuel M. Guadalajara 
Diéguez Y (kN/m3) 135 | 145 | 16.00 | 17.30 
Ts (s) 30m 0.37 
Profundidad (m) | 0-15 | 1.5 -4.5 | 4.5-7.5 | 7.5-60 >60 
Eulogio Parra y 3 Vs (m/s) 95 215 255 570 1000 
11 Guadalajara 
Pablo Casals ; y (kN/m3) 125 14.3 14.50 16.40 | 17.10 
Ts (s) 30m 0.53 
Profundidad (m) | 0-2 2-9 9 - 34 > 34 
do tablo Neruda y a Vs (m/s) 158 220 420 950 
aseo uadalajara 
Jacarandas y (KN/m3) 13.1 14.5 16.10 17.00 
Ts (s) 30m 0.41 
Puno sao pauio a O 
Prol. Av. i S (m/s 
13 ae Guadalajara 
Américas y Sao i y (kN/m3) 125 | 145 15.20 | 17.20 
Paulo Ts (s) 30m 0.16 
Profundidad (m) | 0-2 | 2-6.5 | 6.5-10 | 10-17 > 17 
i Vs (m/s 123 180 250 612 1900 
44 Patria y Eva Zapopan (m/s) 
Briseño Y (kN/m3) 12.9 14.5 14.70 16.50 17.30 
Ts (s) 30m 0.26 
Profundidad (m) | 0-4 4-14 14-46 > 46 
Patria y Paseo Vs (m/s) 180 448 513 1200 
1 Z 
? Royal Country apopan y (kN/m3) 14.50 | 15.40 16.10 17.10 
Ts (s) 30m 0.43 
Profundidad (m) | 0-2 2-5 5-10.5 | 10.5-38 | >38 
Vs (m/s 111 196 607 696 1200 
16 Paseo Valle Real Zapopan (m/s) 
y Federalistas y (kN/m3) 12.5 14.6 16.30 16.40 | 17.10 
Ts (s) 30m 0.33 
Profundidad (m) | 0 - 2.5 | 2.5 -8.5 | 8.5 - 23 23 -72 >72 
iféri Vs (m/s 216 289 445 624 1900 
47 Periférico y Zapopan (m/s) 
Laureles y (kN/m3) 14.60 | 14.70 16.20 16.60 | 17.40 
Ts (s) 30m 0.57 
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Número Ubicación Municipio | Estratos 1 2 3 4 5 6 
ter is A 
ábli s (m/s 
18 Pública, Zapopan 
leal y (kN/m3) 13.10 | 14.70 | 15.50 15.50 | 17.00 
SRT AaS Ts (s) 30m 0.47 








4.1.1 Interpolación de EPU 


Para obtener efectos de sitio, son necesarias las aceleraciones espectrales 
en terreno firme de un punto en específico, en este caso para los puntos de los 
perfiles de suelo (tabla 4). 


Se realizó una interpolación utilizando el método de Ponderación de 
Distancia Inversa (IDW), de Shepard (1968), para cada punto de los 18 mostrados 
en la tabla 4, se determinaron sus Espectros de Peligro Uniforme en terreno firme 
con la ayuda de la malla de interpolación desarrollada en el capítulo 3. Se obtuvieron 
los EPUs para cuatro períodos de retorno, 100 años, representando aceleraciones 
para estructuras tipo C, 475 años, para estructuras tipo B, 975 años para estructuras 
tipo A, y 2475 años para estructuras especiales (estructuras de mayor importancia). 


En la figura 17 se observa la localización del punto 10, “Colomos y Manuel 
M. Diéguez”, respecto a los 24 puntos de la malla de interpolación, y el resultado de 
las aceleraciones espectrales obtenidas por la interpolación entre los puntos ZMG- 
17, ZMG-18, ZMG-12 y ZMG-13 en la figura 18. 





Figura 17. Localización del punto 10, “Colomos y Manuel M. Diéguez” respecto a 
los puntos de la malla de interpolación (Buenrostro, 2017) 
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Figura 18. EPU punto 10 “Colomos y Manuel M. Diéguez”, producto de una 
interpolación (IDW) (Buenrostro, 2017) 


4.1.2 Ajustes de acelerogramas a los Espectros de Peligro Uniforme 


Después de modelar los estratos del suelo, se debe aplicar un movimiento 
o excitación en la base, en este caso en terreno firme. El acelerograma que se usará 
es el registrado en la estación CICEJ Pozo 35m, con aceleraciones máximas de 
6.46, 5.96 y 6.96 cm/s?, para los canales N90W, V y NOOW, respectivamente, 
mostrado en la figura 16. Como se observa las aceleraciones son pequeñas por lo 
que, la señal se debe de escalar al EPU, de cada punto deseado. Para este estudio 
se realizaron dos técnicas de ajustes de espectros y se compararon, para decidir, 
cuál sería el más apropiado. El primer método fue la técnica de Ajuste de Espectros 
(Spectral Match), el cual implica modificar el contenido de las frecuencias de los 
registros sísmicos utilizados (inicial) o registros semilla, hasta que logren empatar 
con el espectro deseado (objetivo) teniendo así un espectro y/o acelerograma 
ajustado. Esta metodología se desarrolló por Al Atik y Abrahamson (2010), la cual 
se implementó en el programa de cómputo RspMatch 2009. Para esta investigación 
se utilizó el apartado de Spectral Matching del programa EZ-FRISKv7.65'M, donde 
se incluye la metodología mencionada en el párrafo anterior. En la figura 19, se 
observa una comparación entre las dos funciones transferencia, obtenidas por los 
dos métodos de escalamiento usados para un período de retorno de 475 años, para 
el perfil de suelo del punto 10 de los 18 antes mencionados, “Colomos y Manuel M. 
Diéguez”. 


Se puede apreciar que los resultados son similares, pero en el ajuste de 
espectros que se hizo por la técnica de ajuste de espectro “spectral matching” (color 
verde), se observa que es ligeramente menor la respuesta, debido a que en esta 
técnica lo que se cuida es el contenido de frecuencias más que las aceleraciones, 
por lo que existe una pequeña pérdida de datos. Se llegó a la conclusión de usar 
para los 18 perfiles de suelo que se tienen de la ZMG, el segundo método de 
escalamiento por amplitud. Usando para todos los casos el acelerograma registrado 
en la estación CICEJ Pozo 35m, en las dos componentes horizontales N-S y E-O. 
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Una vez obtenidas las funciones de transferencia diferentes para cada EPU, 
con diferentes períodos de retorno, se multiplicaron directamente con el EPU (IDW) 
en terreno firme en términos de frecuencias, para así conseguir los espectros de los 
18 puntos afectados con el sitio. 


En la figura 20 se muestran los espectros de sitio transparentes o EPUs con 
efectos de sitio para el punto 10, “Colomos y Manuel M. Diéguez”, con la técnica de 
escalamiento elegida, para períodos de retorno de 100, 475, 975 y 2475 años, en el 
sentido N-S (el mayor registrado). 
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Figura 19. Comparación de las funciones de transferencia para el punto 10, 
“Colomos y Manuel M. Diéguez” 
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Figura 20. Espectros de respuesta elásticos en la superficie para el perfil de 
“Colomos y Manuel M. Diéguez”, registrados en la dirección N-S 
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5 ANÁLISIS COMPARATIVO CON ESPECTROS DE DISEÑO 
5.1 Comparación con la reglamentación actual del sitio estudiado 


Se analizó la reglamentación actual de la ciudad de Guadalajara 
(ROMG,1997), y la del municipio de Zapopan (RCDUMZ, 2003) y se realizó una 
serie de comparaciones de estos reglamentos, contra los Espectros de Diseño en 
terreno firme y en superficie obtenidos en este trabajo, para los puntos de la 
“Catedral de Guadalajara” (Latitud: 20.677, Longitud: -103.347), y “Colomos y 
Manuel M. Diéguez” (Latitud: 20.692, Longitud -103.371), ambos ubicados en el 
municipio de Guadalajara. 


Los EPUs en terreno firme y en superficie (filtrados), para ejemplificarlos 
como un espectro de diseño, se les aplicó una meseta y curvas ascendentes y 
descendentes, también fueron afectados por un factor de sobrerresistencia de 2, ya 
que la comparación se hace con espectros de diseño, que, aunque estos últimos se 
obtengan sin un Q’, eso no nos indica que sean totalmente “transparentes”, ya que 
tienen la sobrerresistencia implícita. La comparativa contra los espectros en 
superficie “filtrados” contra el ROMG,1997 para el punto 1 “Catedral de 
Guadalajara”, se pueden apreciar en la figura 21 donde las caídas de las curvas de 
los espectros propuestos por este análisis son controladas a diferencia de las del 
Reglamento de Guadalajara. 


Se obtuvieron, los espectros de desplazamientos para los espectros de 
diseño propuestos y para los del ROMG, 1997 (ver figura 22), donde son notables 
las diferencias, ya que los espectros del ROMG, 1997 tienden a infinito, y como 
consecuencia dan mayores desplazamientos en las edificaciones. Se controlaron 
las caídas de las curvas con los parámetros “k” y “r”, en los espectros de diseño 
propuestos. Estas comparativas también se analizaron para terreno firme tipo l. 


eoeeeeeeeee Este estudio. Estructuras tipo B. "Tr=475 años" 
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Figura 21. Comparación de espectros de diseño del ROMG, 1997, con los 
espectros de diseño obtenidos en superficie (propuestos) del punto 1, “Catedral de 
Guadalajara”, para suelos tipo Il 
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Figura 22. Comparación de espectros de desplazamientos del ROMG, 1997, con 
los espectros de desplazamientos obtenidos en el punto 1, “Catedral de 
Guadalajara”, para suelos tipo Il 


Se trabajó en la segunda comparación con el Reglamento Orgánico del 
Municipio de Guadalajara, 1997 y el punto 10, “Colomos y Manuel M. Diéguez”, para 
los mismos tipos de estructuras y suelos que en el punto anterior, ya que sus 
características geotécnicas son equiparables con las características de los suelos 
terreno ll que describe el ROMG, 1997. También se realizaron algunas 
comparaciones con los Espectros de Diseño del RCDUMZ, 2003 (Zapopan), para 
estructuras tipo “A3”, “A1 y A2” y “B”, de terreno tipo | y los Espectros de Diseño en 
terreno firme y en superficie para el punto 16, “Paseo Valle Real y Federalistas” 
(Latitud: 20.7209, Longitud: -103.4361), para estructuras tipo “Especiales”, “A” y “B”. 


En las tablas 5 y 6 se pueden notar los PGA de los Espectros de Diseño 
(ED) reducidos en roca y en superficie para el punto 1, “Catedral de Guadalajara” y 
el punto 10, “Colomos y Manuel M. Diéguez”, contra el Reglamento Orgánico 
Municipal de Guadalajara, 1997, para los terrenos | y Il, respectivamente. Donde se 
puede observar que, son muy parecidos los valores en roca o terreno | de ambos 
puntos, pero, es importante tomar como referencia el punto 10, ya que tiene 
amplificaciones mayores, que los de la Catedral, debido a la amplificación de las 
ondas del suelo, y a pesar de ello ambos se catalogarían como un terreno ll en el 
ROMG, 1997. 


Tabla 5. Comparación de PGA, “Catedral de Guadalajara” contra el Reglamento Orgánico 
del Municipio de Guadalajara, 1997 


PGA en roca “Catedral de Guadalajara” contra PGA en superficie, Tr = 475 años “Catedral de 
“ROMG, 1997” terreno | Guadalajara” contra “ROMG, 1997” terreno ll 


RONMG, 1997 ROMG, 1997 
(Terreno I) (Terreno ll) 


0.725 
Tipo A/A1, A2 Tipo A/A1, A2 
Tipo Especiales/A3 0.468 Tipo Especiales/A3 0.468 


PGA(g) 
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Tabla 6. Comparación de PGA, “Colomos y Manuel M. Diéguez” contra el Reglamento 
Orgánico del Municipio de Guadalajara, 1997 


PGA en roca “Colomos y Manuel M. Diéguez” PGA en superficie, Tr = 475 años “Colomos y Manuel M. 
contra “ROMG, 1997” terreno | Diéguez” contra “ROMG, 1997” terreno ll 


ROMG, 1997 ROMG, 1997 
Terreno | Terreno ll 


0.21 
Tipo A/A1, A2 0.26 Tipo A/A1, A2 
Tipo Especiales/A3 | 0.32 0.468 Tipo Especiales/A3 0.468 





También se pueden apreciar en la tabla 7, los PGA del punto 16, “Paseo 
Valle Real y Federalistas” en roca y en superficie contra la normatividad vigente del 
municipio de Zapopan, RCDUMZ, 2003, para terreno I. El punto de Paseo Valle Real 
y Federalistas, es muy ambiguo ya que los espectros en roca y en superficie entran 
en la misma categoría de suelo, según la clasificación de los reglamentos de 
Guadalajara y de Zapopan, por lo cual la comparativa se realiza solo para terreno l. 


Tabla 7. Comparación de PGA, Paseo Valle Real y Federalistas contra el Reglamento de 
Construcciones y Desarrollo Urbano del Municipio de Zapopan, 2003 


PGA en roca “Paseo Valle Real y Federalistas” PGA en superficie, Tr = 475 años “Paseo Valle Real y 
contra “RCDUMZ, 2003” terreno | Federalistas” contra “RCDUMZ, 2003” terreno | 


RCDUMZ, 2003 RCDUMZ, 2003 
PGA(g) Ed (Terreno I) PGA(9) CO (Terreno I) 
loss | 


Tipo A/A1, A2 Tipo A/A1, A2 
Tipo Especiales/A3 0.468 Tipo Especiales/A3 0.468 





5.2 Comparación con PRODISIS del Manual de Obras Civiles de la 
CFE, 2015 


Se obtuvieron EPUs en terreno firme en el programa de cómputo 
PRODISISv4.1, de CFE, 2015, para períodos de retorno de 100 años, 475 años, 
975 años, y 2475 años, para cada uno de los tres sitios anteriormente comparados. 
Los parámetros sísmicos, varían para cada punto y para cada período de retorno, 
por lo que los espectros que se obtengan con este programa de cómputo, son 
directamente comparados, con los EPUs en terreno firme o roca de este estudio. A 
ninguno de los EPUs se les aplicó alguna sobrerresistencia ya que los resultados 
que arroja el PRODISISv4.1, son “transparentes”, y/o “elásticos”. En la tabla 8, se 
muestra ésta comparativa. Donde es posible notar que los EPUs en roca obtenidos 
en PRODISIS, son menores que los conseguidos en este estudio. A pesar de que 
en las recomendaciones del MDOC 2015, se indica que los espectros de respuesta 
en roca para un período de retorno específico, son espectros realizados con un 
enfoque probabilista, la diferencia de las ordenadas espectrales podría darse, en la 
metodología y criterios utilizados en la obtención de estos EPUs. 
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Tabla 8. Comparación de PGA, de EPUSs elásticos contra el PRODISISvA. 1 


PGA en EPUs transparentes PGA en EPUs transparentes PGA en EPUs transparentes 
en roca “Catedral de en roca “Colomos y Manuel en roca “Paseo Valle Real y 
Guadalajara” contra M. Diéguez” contra Federalistas” contra 

“PRODISIS, 2015” “PRODISIS, 2015” “PRODISIS, 2015” 


PRODISI PRODISI PRODISI 
a e [PEE o [o [MEE vu e EE 


Tr =100 


T=475 años Tr=a75 años TAS años 
Tr=975 años Tr=975 años Tr=975 años 


Tr=2475 Tr=2475 Tr=2475 
E E DIKS 


6 CONCLUSIONES 





Se estudió el peligro sísmico para la Zona Metropolitana de Guadalajara, 
ubicada en el estado de Jalisco. Se realizó una microzonificación, mediante una 
malla, con intersecciones a cada 10 km, en donde cada intersección es un punto. 
Para cada punto se obtuvo un Espectro de Peligro Uniforme (EPU), para diferentes 
períodos de retorno y/o probabilidades de excedencia. Se asoció un EPU de un 
período de retorno específico, para el diseño de un tipo de estructura, quedando un 
Tr=100 años como equivalente a un espectro de diseño de estructuras de tipo “C”, 
un Tr=475 años a las estructuras tipo “B”, un Tr=975 años a las estructuras tipo “A” 
y las “Especiales”, equivalentes a un Tr=24'75 años. 


Se realizaron comparaciones para tres puntos (Catedral de Guadalajara, 
Colomos y Manuel M. Diéguez y Paseo Valle Real y Federalistas), de los 18 puntos 
obtenidos con efectos de sitio, con la normatividad vigente del municipio de Zapopan 
y Guadalajara. Para poder hacer una comparación directa se tomaron los espectros 
de los reglamentos (RCDUMZ, 2003 y ROMG, 1997), sin un factor de ductilidad, y 
se compararon con los obtenidos de los EPU en superficie. Se tomaron tres puntos 
en donde la respuesta de aceleración espectral en la superficie fuera diferente. Los 
tres espectros elegidos, se utilizaron como ejemplo para reducirlos por un valor de 
sobrerresistencia de 2 y aplicarle una meseta, la cual cubre las aceleraciones 
máximas del espectro y el período del suelo. De esta forma se pudieron comparar 
los espectros de diseño de los reglamentos con los espectros reducidos. Las 
comparaciones también se realizaron con el programa de cómputo PRODISISv4.1, 
2015. Se utilizaron los EPU en terreno firme elásticos, de los tres puntos, para poder 
ser comparados con los espectros que arroja PRODISISv4.1. 


Se observó que las respuestas espectrales en terreno firme o terreno | para 
los reglamentos de Zapopan y Guadalajara son mayores que los mostrados en éste 
estudio, pero que para un terreno tipo ll, es todo lo contrario, ya que los espectros 
de las normativas no distinguen en aceleración según cambie de terreno, sino, sólo 
amplifica el tamaño de la meseta, dejando una gran incertidumbre en la 
amplificación de los suelos de la ZMG. 


Se puede apreciar que la clasificación de los suelos del RCDUMZ, 2003 y 
ROMG, 1997, no es muy apropiada, ya que pueden entrar varios tipos de suelos 
(con características dinámicas muy diferentes) en un mismo terreno, dejando abierta 
la opción de hacer diferentes interpretaciones del reglamento y corriendo el riesgo 
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de subdiseñar o sobrediseñar las edificaciones de la ciudad. Lo ideal es actualizar 
la clasificación de suelos, con base en estudios geotécnicos y geofísicos recientes, 
y dejar de lado esta clase de incongruencias. 


Con las comparativas realizadas se puede apreciar que los espectros de 
diseño de ambos reglamentos analizados en esta investigación (RCDUMZ, 2003 y 
ROMG, 1997), necesitan una actualización y este estudio ofrece parámetros para 
continuar con la investigación y proponer espectros y parámetros de diseño sísmico 
adecuados para la Zona Metropolitana de Guadalajara. 
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